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O ' Abstract. We compute relativistic corrections to the process e+e J /xj/ + rjc and find that they 

I resolve the discrepancy between theory and experiment. 

(N : 

CS) ■ The disagreement between theory and experiment for the exclusive double- 

. charmonium process e^e^ J/xir + rjc has, for a number of years, been one of 

^ I the largest discrepancies in the standard model. The production cross section times the 

O ' branching fraction into two or more charged tracks o{e^e^ ^ J/\ir + rjc) x 5>2, has 

^ . been measured by the Belle Collaboration to be 25.6 ± 2.8 ± 3.4 fb (Ref. [1]) and by 

! the BABAR Collaboration to be 17.6 ± 2.8+2 i ^ contrast, nonrelativistic 

O I QCD (NRQCD) factorization fs*] calculations at leading order in predict cross 

i!^ ■ sections of 3.78 ± 1.26 fb (Ref. [4]) and 5.5 fb (Ref. [5]). The differences between 

Ph, these calculations arise from different choices of nic, NRQCD matrix elements, and as 

{iLii and from the fact that the calculation of Ref. [|4|] includes QED effects, while that of 

^ ' Ref. riHl] does not. An important recent development is the calculation of the corrections 

■ of next-to-leading order in a^, (Ref. [6]), which yield a K factor of about 1.96. However, 
. !^ i even if one includes this K factor, a significant discrepancy remains. 

^ ■ It is known from the work of Ref. ||4|] that relativistic corrections to the process are 

■ potentially large. The first relativistic correction appears at order v^, where v is the heavy 
quark or antiquark velocity in the quarkonium rest frame, (v^ ~ 0.3 for charmonium.) 
In Ref. [41], the order-v^ corrections are estimated to give a K factor 2.0+j j. The large 
uncertainties arise because of large uncertainties in the relevant nonperturbative NRQCD 
matrix element of relative order v^. If the order- K factor is indeed large, then this 
casts doubt on the convergence of the v expansion of NRQCD. In what follows, we 
address the uncertainty in the order- matrix through a potential-model calculation and 
the convergence of the v expansion through resummation. 
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The NRQCD matrix element of leading-order in v for production (or decay) of the rfc 
is related to the wave function at the origin i//(0) : 

Here i// annihilates a heavy quark, annihilates a heavy antiquark, and V^(x) and V/(p) 
are coordinate-space and momentum- space Schrodinger wave functions, respectively. A 
similar matrix element appears in NRQCD expressions for 7/i/A production and decay. 
The corresponding NRQCD matrix elements of higher order in are given by 

"^^'"^^^^ ^ / (^^'"'^^^^ = ^(0|Z^(-V^)'>|T],). (2) 

We use the following notation for such matrix elements: (p^") = i/a^^") (0)/i//(0) and 
(v^) = (p^)/m^. Previous attempts to determine \j/^^\0) from phenomenology, from 
lattice measurements, and from the Gremm-Kapustin relation [0] have resulted in large 
uncertainties. Even the sign of i//(^)(0) is not known with great confidence from these 
methods. 

Our strategy is to use a potential-model calculation, with the Cornell potential [8], to 
determine \]/^^\0). Details of this calculation can be found in Ref. [0]. If the model 
potential is the exact static QQ potential, then the errors in the potential model are 
of relative order (Ref. [10]). With an appropriate choice of parameters, the Cornell 
potential provides a good fit to lattice data for the static QQ potential. 

The matrix element i/a^^) (0) contains a linear UV divergence and must be regulated. 
Because existing calculations in NRQCD of order a^. and higher make use of dimen- 
sional regularization, we ultimately want to obtain a dimensionally regulated matrix 
element. We introduce two methods to achieve this. 

Method 1 requires knowledge only of the wave function. In this method, we first 
regulate using a simple, analytic momentum-space hard cutoff A^/ {p^ + A^). Then, we 
calculate the difference between (0) in hard-cutoff and dimensional regularization, 
which we call Av/(2)(o). Since Av/(2)(o) gives the difference between UV regulators, it 
is dominated by large momenta and can be computed in perturbation theory. We subtract 
Ai//(^)(0) from the hard-cutoff result. In order to obtain the dimensionally regulated 
result, we must extrapolate this subtracted expression to A = oo because our computation 
of Ai//(^) (0) does not include contributions that are suppressed as 1/A. 

Method 2 is applicable only to potential models. It requires knowledge of the binding 
energy and the potential, in addition to the wave function. In this approach, we use 
the Bethe-Salpeter equation to expose an explicit loop from the wave function in the 
expression for the matrix element. Then, we regulate that loop dimensionally. The result 
is equal to that which would be obtained from the Gremm-Kapustin relation, but for the 
binding energy of the potential model. 

The results from Method 1 and Method 2 agree well numerically and yield i/a^^) (0) = 
0.1 18 ±0.024 ±0.035 GeV^/^ ^^d (p2) = 0.50 ±0.09 ±0.15 GeV^, which imply for 
mc = \A GeV that (v^) ^ 0.25 ± 0.05 ± 0.08. This last result is in good agreement with 



expectations from the NRQCD v-scaling rules. In all of these results, the first error bar 
reflects the uncertainty in the input potential-model parameters and the wave function at 
the origin, and the second error bar reflects relative-order- corrections that have been 
neglected. This is the first determination of l/A*^^) (0) with small enough uncertainties to 
be useful phenomenologically. 

We can extend Method 2 to matrix elements of higher order in v. By using the equation 
of motion, dimensional regularization, and the scalelessness of the individual terms in 
the Cornell potential, we obtain the simple relation 

= (p2)". (3) 

This relation allows one to resum a class of the relativistic corrections to 5-wave quarko- 
nium decay and production amplitudes to all orders in v. 

We now apply these results to the relativistic corrections to a[e+e^ ^ 7/1//+ T]c]. 
These corrections arise in two ways. First, they appear directly in the process e+e^ 
J/\j/ + ric itself. Second, they enter indirectly through 1 1//(0) p, which appears as a factor 
in the contribution of leading order in v^. The quantity 1 1//(0) p is determined from the 
experimental value for the width for J/xjr ^ e+e^ and the theoretical expression for that 
process, which is affected by relativistic corrections. 

Our preliminary results, including the effects of QED, as well as QCD, contributions 
are as follows. The resummed relativistic corrections to o[e^e^ J/xir + rjc] itself 
yield a K factor 1.34, while the resummed relativistic corrections to |i//(0)p yield a 
K factor 1.32. In both cases, the resummation of contributions beyond order has 
only about a 10% effect on the K factors, which indicates that v expansion of NRQCD 
converges well for these processes. Combining the resummed relativistic-correction 
K factors with the relative-order- K factor 1.96, we obtain a complete K factor 
4.15. Applying this to the calculation of Ref. ^ and using the most recent value for 
r[7/v/^e+e"] (Ref. [1 1]), we obtain aprediction C7[e+e" ^y/i/A-t-r]^ = 17.5±5.7fb. 
The quoted uncertainty reflects only the uncertainties in the values of and (p^) . Other 
theoretical uncertainties are large and need to be quantified. Nevertheless, our prediction 
is in agreement with the Belle and BABAR results, and it seems that the inclusion of 
relative-order- a^, and relativistic corrections resolves the discrepancy between theory 
and experiment at the present level of precision. 
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